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This  paper  presents  the  oxygen  stoichiometry  control  problem  of  proton  exchange  membrane  (PEM) 
fuel  cells  and  introduces  a  solution  through  an  optimal  control  methodology.  Based  on  the  study  of  a 
non-linear  dynamical  model  of  a  laboratory  PEM  fuel  cell  system  and  its  associated  components  (air 
compressor,  humidifiers,  line  heaters,  valves,  etc.),  a  control  strategy  for  the  oxygen  stoichiometry  reg¬ 
ulation  in  the  cathode  line  is  designed  and  tested.  From  a  linearised  model  of  the  system,  an  LQR/LQG 
controller  is  designed  to  give  a  solution  to  the  stated  control  problem.  Experimental  results  show  the 
effectiveness  of  the  proposed  controllers  design. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

From  the  automatic  control  point  of  view,  a  fuel-cell-based  (FCB) 
system  is  a  nonlinear  dynamic  plant,  with  multiple  inputs,  multiple 
outputs,  variables  strongly  coupled,  model  uncertainty  and  inci¬ 
dence  of  external  disturbances.  In  this  context,  a  topic  that  deserves 
special  attention  is  the  oxygen  stoichiometry  control  [1-3].  If  the 
oxygen  flow  at  the  cathode  of  a  fuel  cell  is  too  low,  it  produces  hot 
spots  on  the  polymeric  membrane,  decrementing  the  cell  power 
due  to  the  lack  of  reactant  in  the  triple  contact  areas.  This  phe¬ 
nomenon  is  known  as  cathode  starvation  and  can  lead  to  irreversible 
damages  in  the  membranes  [1,4].  On  the  other  hand,  an  exces¬ 
sively  high  oxygen  flow  could  cause  a  fast  drying  phenomenon 
of  the  polymeric  membranes,  increasing  its  ionic  resistance.  Apart 
from  that,  higher  air  flows  mean  more  compressor  power  consump¬ 
tion,  degrading  the  overall  system  energy  conversion  efficiency. 
Therefore,  in  order  to  avoid  irreversible  damages  in  the  polymeric 
membranes  and  to  obtain  an  efficient  and  reliable  power  response, 
it  is  necessary  to  design  a  control  system  capable  of  regulating 
the  oxygen  stoichiometry  and  thereby  preserve  proper  operation 
conditions  in  the  system. 

According  to  preliminary  works  presented  in  [1,5],  the  linear 
quadratic  regulator  (LQR)  approach  represents  a  promising  option 
to  start  solving  the  stoichiometry  control  problem  in  PEM  FCB 
systems  with  air  compressors.  In  these  first  approaches,  some  ini¬ 
tial  assumptions  such  as  ideal  static  humidifiers  were  taken  into 
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account.  On  the  other  hand,  the  LQR  controller  law  is  based  on  the 
availability  of  the  complete  state  vector,  which  is  not  fully  measur¬ 
able  in  all  cases.  To  overcome  this  inconvenient,  an  observer  design 
based  on  three  measurements  (stack  voltage,  compressor  air  flow 
and  manifold  pressure)  was  developed  to  estimate  the  state  vector. 
Finally,  all  the  control  tests  were  performed  through  simulations 
with  the  plant  non-linear  model. 

The  current  control  strategy,  presented  in  this  paper,  resumes 
the  idea  firstly  proposed  by  Pukrushpan  et  al.  using  an  LQR  design 
for  air  regulation  and  system  optimisation.  Nevertheless,  the  non¬ 
linear  model  of  the  plant  was  completely  validated  in  the  test  bench 
[6],  the  air  humidifier  was  modelled  as  a  dynamical  system  to 
include  its  slow  response  and  a  reduced-order  Kalman  observer 
based  on  single  measurement  was  designed  to  estimate  the  two 
states  needed  for  the  control  algorithm.  Furthermore,  the  controller 
strategy  was  implemented  and  the  parameters  retuned  in  the  fuel 
cells  test  bench  [7]. 

Therefore,  the  main  contributions  of  this  paper  lead  to  the  appli¬ 
cation  of  LQR  designs  in  the  compressor  air  flow  taking  into  account 
several  scenarios  related  to  variations  in  both  the  oxygen  stoi¬ 
chiometry  and  the  stack  current.  Although  those  designs  were 
performed  employing  the  linearisation  of  the  validated  mathemat¬ 
ical  model  in  [6],  different  controller  set-ups  have  been  stated  and 
implemented  over  a  real  FCB  system.  Experimental  results  demon¬ 
strate  the  effectiveness  of  the  control  law  for  all  the  considered 
scenarios.  Additionally,  practical  issues  related  to  the  implemen¬ 
tation  process  and  the  fuel  cell  test  bench  are  presented  and 
discussed. 

The  reminder  of  this  paper  is  organised  as  follows:  in  Section  2, 
a  detailed  description  of  the  real  FCB  system  is  presented  as  well 
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Nomenclature 

System 

TM hum,ca 

mass  of  air  in  the  cathode  humidifier  [kg] 

1TIv2  > an 

mass  of  vapour  in  the  stack  anode  [kg] 

PO.V2 ,  ca 

mass  of  vapour  in  the  stack  cathode  [kg] 

2,  an 

mass  of  hydrogen  in  the  stack  anode  [kg] 

^0  2,ca 

mass  of  oxygen  in  the  stack  cathode  [kg] 

2,ca 

mass  of  nitrogen  in  the  stack  cathode  [kg] 

n 

number  of  cells  in  the  stack 

1st 

stack  current  [A] 

F 

Faraday  constant  [C  mol-1  ] 

Mq2 

oxygen  molar  mass  [kg  mol-1  ] 

Vcp 

compressor  DC  voltage  [V] 

Pnet 

output  net  power  [W] 

wcp 

compressor  air  flow  [kg  s-1  ] 

WCp,ref 

desired  compressor  air  flow  [kg  s-1  ] 

WH2 

cathode  hydrogen  flow  [kg  s-1  ] 

^02 

cathode  oxygen  stoichiometry 

■^0  2,des 

desired  cathode  oxygen  stoichiometry 

0)cp 

compressor  shaft  angular  speed  [rad  s-1  ] 

X02 

molar  fraction  of  oxygen  in  the  air 

&  atm 

relative  humidity  of  the  ambient  air 

Controller 

f 

drift  vector  field 

g 

control  vector  field 

u 

control  action  of  the  non-linear  system 

V 

control  action  of  the  linearised  system 

Vp 

noise  vector  of  the  linearised  system 

X 

state  vector  of  the  non-linear  system 

y 

output  vector  of  the  linearised  system 

z 

state  vector  of  the  linearised  system 

z 

observed  state  vector  of  the  linearised  system 

A 

linearised  system  matrix 

B 

linearised  control  matrix 

C 

linearised  output  matrix 

J 

LQR  cost  function 

I<r 

Kalman  estimation  gain 

I<r 

state  feedback  gain 

Q 

LQR  design  matrix 

R 

LQR  design  matrix 

wp 

linearised  perturbation  matrix 

Table  1 

Description  of  the  laboratory  elements  depicted  in  Fig.  2. 


Element  number 

Description 

1 

Electronic  load 

2 

Air  compressor 

3 

Stack  temperature  control  unit 

4 

Cathode  back  pressure  valve 

5 

Air/oxygen  humidifier 

6 

Hydrogen  humidifier 

7 

PEM  fuel  cell  stack 

8 

Real  time  computer 

9 

Stack  heater 

10 

Temperature  sensor 

11 

Current  collector 

12 

Bipolar  plates 

assemblies  (MEAs),  with  a  catalyst  loading  of  1  mg  cm-2  of  plat¬ 
inum,  50  cm2  of  active  area,  50  W  of  nominal  power  and  100W 
peak  power.  Cellkraft®  membrane  exchange  humidifiers  are  used 
to  maintain  proper  humidity  conditions  inside  the  cells,  which  is 
crucial  to  ensure  the  optimal  operation  of  PEM  membranes.  The 
air  compressor  consists  in  a  12  V  DC  oil-free  diaphragm  vacuum 
pump.  The  line  heaters  and  stack  temperatures  are  controlled  by 
using  a  power  station  via  decentralized  PID  controllers,  allowing 
independent  gas  conditions  (humidity  and  temperature)  inside  the 
stack. 

2.2.  Nonlinear  model  of  the  system 

Adopting  a  nonlinear  control  design  approach,  all  the  subsys¬ 
tems  have  been  modelled  and  analyzed  individually  in  order  to 
obtain  a  nonlinear  state-space  model  [6].  After  regrouping  the 
resultant  differential  equations,  in  conjunction  with  the  static 
relationships  between  variables,  the  following  description  of 
seventh-order  nonlinear  control  affine  system  was  obtained: 

^p=/(x(t))  +  g«t))u(t),  (1) 

where  x  e  R7+  corresponds  to  the  vector  of  the  states  variables  and 
ueR+  denotes  de  control  variable.  The  description  of  the  state  vari¬ 
ables  and  the  control  action  is  outlined  as  follows: 

^1  —  0)cp,  X2  =  1Tlhum,cai  ^3  —  ^Zo2,  ca-> 

X4  =  nTN2?ca>  *5  —  TTlv2,cai  *6  —  mH2,an? 

X-j  =  mV2,ani  tl  =  Vcp- 


as  the  system  model  considered  for  controller  designs.  Section  3 
briefly  explains  the  LQR/LQG  strategy  and  the  considered  control 
objectives.  In  Section  4  the  most  relevant  results  obtained  from 
the  implementation  of  the  controller  designs  are  presented  and 
discussed.  The  main  conclusions  and  some  lines  of  further  research 
are  exposed  in  Section  5. 

2.  PEM  fuel  cell  test  bench  and  control  model 

2.1.  System  under  study 

The  laboratory  system  under  consideration  is  mainly  composed 
by  a  fuel  cell  stack,  an  air  compressor,  a  hydrogen  storage  tank, 
gases  manifolds,  humidifiers  and  line  heaters.  The  schematic  dia¬ 
gram  in  Fig.  1  depicts  the  PEM  fuel  cell  laboratory  arrangement, 
Fig.  2  shows  the  real  test  bench  elements  and  Table  1  presents  their 
description. 

Some  of  the  technical  specifications  of  the  laboratory  fuel  cell 
system  under  study  are  summarised  as  follows:  the  stack  is  an 
ElectroChem®  7-cell  stack  with  Nation  115®  membrane  electrodes 


2.3.  Model  linearisation 


In  order  to  design  LQG  controllers  for  the  system  under  study,  it 
is  necessary  to  preform  a  model  linearisation  in  a  neighbourhood 
of  a  given  working  point.  The  linearisation  procedure  consists  in 
expanding  the  right  side  of  (1 )  in  Taylor  series  around  a  stationary 
operating  point  and  consider  its  first  derivative  (Jacobian  matrix). 
The  choice  of  the  working  point  was  based  on  practical  consider¬ 
ations  and  in  a  way  to  get  an  average  point  within  the  operation 
range  of  the  stack.  The  Taylor  development  of  (1 )  around  x(0)  =  x0 
at  the  first  order  in  z(t)  =  x(t)  -x(0)  is  given  by 


dzjt) 

dt 


mm) 

am 


z  + 

x(t)=x0 


dg(m) 

dm 


(u(t)-uo) 

X(t)=X0 


Az(t)  +  Bv(t), 


(2) 


where  u0  =  u( 0)  and  v(t )  =  u(t )  -  u0.  The  working  point  selected  in 
this  paper  was  defined  considering  a  compressor  voltage  ( Vcp  =  u0) 
of  6  V,  a  load  current  (Jst)  2  A,  anode  and  cathode  pressure  of  1  bar 
and  a  hydrogen  flow  (  Wh2  )  of  2  slpm.  This  fact  implies  the  following 
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C  ^  Measured  variable 
f  **  Control  action 


Load 


©1 


:'Vcp';. 
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Air  Compressor 
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dp 

o2+  n2+  h2o 

Hum.  -  Cathode  LH 

1  ^ - - ^ 
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■CSD 
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Fuel 
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Stack 


IK 

Back  pressure 
regulators 
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Fig.  1.  Schematic  diagram  of  the  experimental  test  station. 


Fig.  2.  Fuel  cell  test  bench  at  IRI  (CSIC-UPC). 
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values  for  the  elements  of  the  initial  condition  vector: 


*o  = 


207.3136 
2.558  x  10"4 
4.1  x  10“4 
1.415  x  10“3 
1.862  x  10"4 
1.27  x  10“4 
1.844  x  10“4 


3.  Controller  design 

3.1.  Control  objectives 


following  cost  function  is  minimized: 

J  =  lim  l  [  [z(t)T  Qz{t)  +  v(t)TRv(t)]dt,  (4) 

r^o°T  Jo 

where  the  tuning  matrices  Qand  R  are  properly  selected,  fulfilling 
the  design  conditions  Q=Qr>0  and  R  =  RT>  0.  These  matrices  are 
known  as  the  input  I  state  penalization  matrices ,  respectively,  which 
are  in  charge  of  the  prioritization  of  either  the  state  tracking  or  the 
control  energy.  The  solution  of  the  LQR  problem  consists  in  finding 
a  stabilizing  gain  I<r  for  the  state  feedback  law 

v(t)±-Krz(t),  (5) 

where  Kr  acts  as  the  gain  matrix  of  the  controller  in  the  closed-loop 
scheme  (for  further  details  please  refer  to  [9]). 


In  the  world  of  analysis  and  control  of  FCB  systems,  there  are 
diverse  control  goals  related  to  different  operational  aspects  of  the 
system,  such  as  the  maximisation  of  its  effective  life,  the  optimi¬ 
sation  of  its  electric  benefits,  the  minimization  of  its  maintenance. 
The  approach  proposed  in  this  paper  is  focused  on  maximizing  the 
output  net  power  of  the  complete  energy  conversion  system  (Pnet), 
i.e.,  the  difference  of  the  electric  power  of  the  fuel  cell  and  the  air 
compressor  consumption. 

As  proposed  in  previous  works  [1,2,8],  this  power  can  be  maxi¬ 
mized  through  the  proper  regulation  of  the  oxygen  stoichiometry 
at  the  fuel  cell  stack.  This  latter  objective  is,  in  turn,  achieved  by 
controlling  the  mass  inflow  of  oxygen  to  the  fuel  cell  cathode. 
This  fact  allows  to  maintain  at  any  time  an  oxidizer  flow  that  can 
be  efficiently  reduced  close  to  the  catalytic  areas  of  the  MEA  and 
hence  satisfying  the  power  requirements  of  the  system  load,  taking 
always  into  account  the  parasitic  consumption  of  the  air  compres¬ 
sor.  However,  several  issues  are  present  in  the  experimental  set-up 
from  the  control  and  instrumentation  point  of  view.  Those  issues 
can  be,  among  others,  the  inaccessibility  of  certain  variables  and 
measurement  noise  characteristic  of  the  sensors. 

In  this  context,  the  design  and  employment  of  control  schemes 
based  on  Linear  Quadratic  Gaussian  (LQG)  laws  are  proposed  as 
promising  options  in  this  framework,  since  these  strategies  col¬ 
lect  simultaneously  the  improvement  of  the  control  objectives  and 
the  robust  estimation  of  the  system  states.  Therefore,  this  section 
proposes  and  discusses  an  LQG  design  oriented  to  the  net  power 
maximization  of  a  FCB  system. 


3.2.  LQR  problem 

The  proposed  approach  regards  to  the  design  of  a  LQR  and  its 
variant  LQG  due  to  the  consideration  of  noises  and  disturbance  of 
Gaussian  spectral  distribution  affecting  the  closed  loop.  The  main 
assumptions  of  this  control  strategy  are  the  statements  of  a  time- 
invariant  linear  system  and  stochastic  noise  models  with  known 
statistic  properties.  Thus,  the  continuous-time  linear  model  of  the 
system  (1),  obtained  through  a  standard  linearisation  procedure 
(see  Section  2.3),  can  be  written  as 

^p-  =  Az(t)  +  By(t)  +  Wplst(t),  (3a) 

y(t)  =  Cz{t)  +  vp,  (3b) 

where  z  e  RUs  is  the  state  vector,  v  e  Rm  denotes  the  vector  of  manip¬ 
ulated  variables,  y  e  Rq  is  the  vector  of  the  system  outputs,  Ist  e  Rp  is 
the  disturbance  vector,  and  vp  e  Rq  corresponds  to  a  vector  con¬ 
taining  the  considered  measurement  noises.  A,  B,  Wp  and  C  are 
the  system  matrices  of  suitable  dimensions.  The  control  strategy 
aims  to  compute  the  set  of  optimal  control  actions  v(t)  such  as  the 


3.3.  State  estimation 

In  order  to  compute  matrix  Kr ,  it  is  necessary  to  have  avail¬ 
able  the  states  Zi  and  z2  (corresponding  to  the  variables  cocp  and 
mhum,ca  of  the  linearised  model),  feature  that  is  not  always  possible. 
Therefore,  it  is  necessary  to  design  a  reduced-order  state  observer 
from  the  measurement  of  the  output  variable  Wcp.  This  auxiliary 
dynamic  system  estimates  the  trajectories  of  these  two  desired 
variables  of  the  fuel  cell  system  with  an  error  and  convergence 
rate  given  by  the  design  parameter  I<e,  which  corresponds  with  the 
estimation  gain. 

Thus,  the  control  law  (5)  is  computed  as  v(t)  =  - I<rz(t ),  where 
z(t)  is  the  reduced  state  estimate  vector  from  the  observer.  Notice 
that  the  controller  and  the  observer  designs  can  be  done  inde¬ 
pendently  (separation  principle).  On  the  other  hand,  in  order  to 
improve  the  closed-loop  performance,  and  specifically  to  eliminate 
the  steady-state  regulation  error,  an  integrator  is  usually  added. 
Hence,  the  design  of  the  LQR  control  law  implies  computing  a  gain 
Kr  considering  the  extended  dynamic  system  [9]. 

4.  Results  and  discussion 

As  formerly  introduced,  the  pursued  control  objective  of  this 
regulation  strategy  is  the  optimization  of  the  energy  conversion  of 
the  FCB  system,  maximizing  the  net  power  generated  by  the  system 
{Pnet)  under  different  load  conditions.  In  Fig.  3,  it  can  be  appreciated 
that  accomplishing  such  optimal  working  condition  is  equivalent 
to  maintain  the  cathode  line  oxygen  stoichiometry  (A,o2 )  close  to  an 
optimal  value.  This  optimal  value  can  be  determined  from  a  thor¬ 
ough  off-line  analysis  of  the  open-loop  system,  considering  changes 
in  the  current  drained  from  the  stack  and  a  wide  set  of  stoichiom¬ 
etry  values  (Fig.  3).  In  an  air-flow  control  scheme,  the  expression 
of  the  compressor  air  reference  can  be  readily  obtained  from  the 


Fig.  3.  System  energy  conversion  efficiency  for  different  load  conditions  ( Pnet  vs. 
A-o2). 
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Control 


Fig.  4.  Schematic  diagram  of  the  controller  structure. 


desired  value  of  the  oxygen  stoichiometry  (A,02  des).  Given  that  the 
molar  fraction  of  oxygen  in  the  air  ( xo2 )  is  a  known  parameter 
and  taking  into  account  a  constant  relative  humidity  of  the  ambi¬ 
ent  air  (£2 atm),  the  desired  mass  flow  of  air  (Wcp  re^)  can  be  directly 
computed  from  the  following  equation: 

Wcp>re/  =  (1  +  i2a[m)-Lx02tdesM02^,  (6) 

where  n  is  the  number  of  cells  in  the  stack,  Mq2  the  oxygen  molar 
mass  and  F  de  Faraday  constant. 


4.1.  Controller  set-up 

The  selection  of  the  matrices  pair  Q-R  for  the  final  controller  was 
done  following  an  iterative  process,  starting  from  the  simulations 
design  and  then  adjusting  the  parameters  in  order  to  fit  them  to 
experimental  constraints.  The  results  of  one  of  the  best  selections, 
from  the  performance  point  of  view,  are  presented  in  the  next  sub¬ 
section.  As  discussed  in  Section  3,  the  system  state  vector  has  been 
expanded  with  an  integrator  in  order  to  eliminate  the  steady-state 
error  of  the  air  flow  regulation.  Besides,  a  classical  anti-windup 
structure  has  been  added  to  prevent  the  winding  of  the  integrator 
when  saturation  occurs  in  the  compressor  actuator.  In  Fig.  4,  the 
implemented  controller  structure  is  depicted. 


4.2.  Experimental  results 

The  objective  of  this  subsection  is  to  demonstrate  the  perfor¬ 
mance  of  the  proposed  LQG  control  strategy  in  the  real  operating 
conditions  and  point  out  some  considerations  about  practical 
issues.  The  first  series  of  tests  was  performed  considering  a  fixed 
oxygen  stoichiometry  and  changes  in  the  stack  current  (Fig.  5).  It 
is  interesting  to  observe  that,  during  the  first  eight  seconds  of  the 
test,  the  control  action  was  saturated  in  its  lowest  value,  but  when 
stack  current  changes  to  2.5  A  the  controller  recovers  its  regular 
behaviour  by  means  of  the  anti-windup  algorithm.  Next,  the  oxy¬ 
gen  stoichiometry  can  rapidly  start  tracking  its  reference.  During 
the  rest  of  the  test,  it  can  be  noticed  the  regulation  behaviour  of 
the  loop  at  different  stack  currents  and  the  dynamic  disturbance 
rejection  ability  of  the  system.  Figs.  5  and  6  exhibit  the  stack  cur¬ 
rent  variation  ( Jst),  the  oxygen  stoichiometry  (\<j2 )  and  its  reference 
(^o2,re/)>  the  motor  voltage  (Vcp),  the  compressor  air  flow  ( Wcp)  and 
the  internal  reference  calculated  by  the  controller  ( Wcpref). 

Afterwards,  a  different  set  of  tests  was  performed  in  the  PEM  fuel 
cell  test  station,  considering  a  constant  load  current  and  oxygen  sto¬ 
ichiometry  reference  variations.  In  Fig.  6,  the  control  performance 
is  shown  when  the  internal  variable  is  regulated  in  a  wide  range 
of  values. 

The  theoretical  analysis  and  the  experimental  results  show 
that  the  adopted  control  strategy  stabilises  the  system  around  the 
working  point  and  endows  an  improved  dynamic  behaviour  to  the 
closed  loop.  An  important  practical  issue  that  needs  to  be  carefully 
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Fig.  5.  Xq2  regulation  for  variable  load  conditions. 
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considered  is  the  integrator  gain  tuning.  Indeed,  a  choice  of  a  high 
value  ensures  a  fast  response  while  causing  permanent  oscillations 
around  the  steady-state  operation  of  the  compressor.  On  the 
other  hand,  a  low  value  ensures  smooth  response  in  steady-state 
operation,  but  presents  degradation  in  the  transient  response.  This 
compromise  can  be  noticed  experimentally. 

5.  Conclusions 

In  this  paper,  a  linear  control  approach  that  solves  an  impor¬ 
tant  control  problem  of  PEM  FCB  generation  systems  is  presented. 
The  control  law  is  based  on  the  proper  design  of  an  LQR/LQG 
strategy.  Such  controller  allows  maximizing  the  energy  conversion 
efficiency  of  the  overall  system  and  avoids  the  oxygen  starvation 
problem.  The  designed  control  strategy  efficiently  regulates  the 
stack  oxygen  stoichiometry  and  requires  reduced  information  to 
operate  since  its  law  only  depends  on  two  easily  measured  vari¬ 
ables  (Wcp  and  1st).  The  control  design  is  based  on  a  linear  model 
of  the  plant,  but  its  analysis  and  validation  is  performed  through 
experimental  tests  for  different  load  conditions.  The  experimental 
results  show  the  effectiveness  and  feasibility  of  the  proposed  con¬ 
trol  approach  given  real  scenarios  of  both  oxygen  stoichiometry 
and  stack  current  variations. 

Currently,  further  research  is  being  performed  by  the  research 
group  in  the  fuel  cells  control  field,  tackling  nonlinear  control 
solutions  (e.g.  higher  order  sliding  modes  and  model  predictive 
control),  different  control  problems  (e.g.  membrane  water  trans¬ 
port  and  temperature  control)  and  other  plant  configurations  (e.g. 
open  cathode  fuel  cells,  self-humidified  systems  and  hybrid  power 
generation  plants). 
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